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ABSTRACT 

Despite the extensive study of lithium depletion during pre-main-sequence con- 
traction, studies of individual stars show discrepancies between ages determined 



from the HR diagrain and ages determined from lithium depletion (jSong et al 



2002 



White &: Hillenbrand! l2005l ) indicating open questions in the pre-main-sequence evolu- 
tionary models. To further test these models, we present high resolution spectra for 
members of the /3 Pictoris Moving Group (BPMG), which is young and nearby. We 
measure equivalent widths of the 6707. SA Li I line in these stars and use them to de- 
termine lithium abundances. We combine the lithium abundance with the predictions 
of pre-main-sequence evolutionary models in order to calculate a lithium depletion age 
for each star. We compare this age to the age predicted by the HR diagram of the same 
model. We find that the evolutionary models under-predict the amo unt of lithium de 



pletio n for the BPMG given its nominal HR diagram age of ~12 Myr (jZuckerman et al 



2001 



, particularly for the mid-M stars, which have no observable Li I line. This results 
in systematically older ages calculated from lithium depletion isochrones than from the 
HR diagram. We suggest that this discrepancy may be related to the discrepancy be- 
tween measured M-dwarf radii and the smaller radii predicted by evolutionary models. 

Subject headings: stars: abundances — stars: pre-main-sequence — stars: late-type 
— stars: evolution — open clusters and associations: individual (/3 Pictoris Moving 
Group) — stars: activity 
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Introduction 



Lithium depletion during pre~main-sequence (PMS) contraction has been extensively studied. 
Studies of open c lusters have shown that lithium depletion is a strong function of both age an d 



stellar mass (e.g. ISoderblom et al.lll993bl : iBarrado v Navascues et al. 



2004 



Mentuch et al.ll2008l l. 



These studies clearly show that the presence of a strong Li I line at 6707. 8A is an indicator of youth 
in late- type stars. The regularity of lithium depletion in op en clusters implies that the lithium line 



can act as a mass-dependent clock ([Jeffries &: Navloiil200ll ). By combining a measurement of the 



lithium abundance of a star with theoretical evolutionary models, one might be able to derive an 
age of a PMS star whose distance is unknown, e.g. a young field star. In order to do this, we need 
to understand PMS lithium depletion very well. 

Most of the work comparing lithium depletion to HR diagram ages has been done in open 
clusters, focusing on comparing the cluster age derived from the lithium depletion boundary 
(ldbA in t he lowest-mass stars with the HR diagram age found from fitting the upper main se- 
quence (e.g. lBarra,do v Navascues. Stauffer. &: Pattenlll999l : IStauffer et al.lll999l : iBurke et al.ll2004l : 



Ijeffries &: Oliveira 



2005 



Manzi et al 



200 



3). Most studies find that the LDB age tends to be older 



than the age calculat ed from upper-main sequence fitting, thoug h the discrepancy may be smaller 
for younger clusters (jjeffries &: Oliveiral l2005l : iManzi et al.l l2008l ) . These results are used to argue 
for the inclusion of convective overshooting in evolutionary models, since it increases the ages de- 
rived from upper-main sequence fitting without affecting the LDB age. These older upper-main 
sequence fitting ages tend to be in better agreement with LDB ages (jBurke et al.ll2004l ). 



A few studies have been done of individual stars, showing that ages derived from comparison 
of lithium abundances with mode l s of l ithium depletion are persistently older than ages derived 
from the HR diagram. ISong et al.l (|2002l ) find that the age of the binary HIP 112312 derived from 
its lithium depletion is > 20 Myr. They find that this system is likely a member of the /3 Pictoris 
Moving Group (BPMG), whose proposed age d erived from pr e-main-seq uence isochr on es on the 
HR diagram is ~ 12 Myr (jZuckerman et al.ll200ll ). Likewise, IWhite k. Hillenbrand! (j2005l ) find that 
the age of St 34 derived from its position on the HR diagram is 8 =b 3 Myr whereas its lithium 
depletion implies an age of > 35 Myr. These results indicate some open questions about lithium 
depletion in the evolutionary models. 

In order to avoid confusing lithium depletion effects with temperature- or mass-dependent 
systematic trends in the models, we would like to compare HR diagram ages with lithium depletion 
ages for individual stars as was done for HIP 112312 and St 34, rather than comparing the ages of 
more massive stars determined from one technique with the ages of less massive stars determined 
from another. To do this, we need a coeval group of late-type stars so we can study lithium depletion 



^For stars that are fully convective before they reach the main sequence [M < O.SA/©), at a given age there is a 
sharp bound ary between stars that have fully deple ted their lithium and slightly less massive stars with no evidence 
of depletion (jSasri et al.lll99d : iBildsten et al.lll997l ). 
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over a range of masses. We want the stars to be of intermediate age, ~ 5-80 Myr; much younger, 
and very little lithium depletion has occurred; much older, and they are on the main sequence, 
where the HR diagram age is degenerate. We also would like them to be nearby, so we can measure 
accurate distances and luminosities, and bright, so we can get high S/N spectra. 

In this paper, we examine lithium depletion in some of the later-type members of the BPMG, 
whi ch is ideal for studying lithium depletion because it is young and nearby (~ 12 Myr, 10-50 
pc; 



Zuckerman et al. 2001 



). We compare the lithium depletion of our sample to the predictions 
of different theoretical PMS models. We also calculate ages from these models for individual 
stars and compare the HR diagram age of each star with the age implied by its lithium depletion 
to see if each model is internally consistent in the two ages determined for a given star. We 
begin by presenting our spectra and lithium equivalent widths in ^ We then determine effective 
temperatures, luminosities, and lithium abundances for the stars in our sample in ^ In SjU we 
compare these data with the HR diagram and lithium depletion isochrones of theoretical models. 
We find that the models systematically under-predict the observed lithium depletion for M stars (or 
equivalently, over-predict their lithium depletion ages compared to HR diagram ages). We suggest 
that the observed tendency for models of main-sequence M stars to under-predict stellar radii may 
hold the key to resolving the discrepancy between the observed and predicted lithium depletion in 
late-type stars. We conclude in S}6l 



2. Observations and Data 

We took spectra o f ten of the late K and M stars or binary systems listed in Table 1 of 



Zuckerman et al.l (120011 ). with separate spectra for each of the components of the binary stars GJ 
799 and BD — 17°6128; the spectra are shown in Figures [1] and [2j We also took a spectrum of a 
somewhat earlier-type member of the BPMG, the spectroscopic binary HD 155555 A/B. The spectra 
were taken with the 4-meter Blanco telescope at CTIO on six consecutive nights between 2001 Oct. 
30 and 2001 Nov. 5, using the echelle spectrograph. The spectrograph covered a wavelength range 
from 4800 A to 8400 A with a measured spectral resolving power of i? ~ 40, 000. The spectra were 
reduced in according to standard procedures. We flux calibrated the spectra with respect 

to two flux standards, HR 9087 and HR 1544, using standard, sensfunc, and calibrate in IRAF. We 
visually examined the flux calibrated spectra to ensure that the calibration was consistent between 
orders and between the two calibration standards. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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2.1. Spectral Type 



To determine spectral types for our sample, we compared the stre ngth of the Ti05 band 
to the strength of the nearby continuum following the method given in iReid et al.l ()1995l ). We 
calculated the ra tio of the flux from 7126A-7135A to the flux from 7042A-7046A. Using the linear 
fit relation from iReid et alJ (119951 ). we t ook the rat i o of these fluxes to determine the spectral 
types of our sample stars. According to iReid et al.l (jl995l ). this method gives an uncertainty of 
±0.5 spectral subtype. Note that this relation is not well-defined for spectral types much earlier 
than MO because the TiO bandheads are very weak or disappear entirely. Thus, the spectral types 
measured for HIP 29964, BD -17°6128 A, and CD -64° 1208 are more uncertain, although probably 
correct in a relative sense. Additionally, since the relation gives a numeric value for the spectral 
type, with representing MO, it is not clear how to assign that number to the Morgan-Keenan 
system for values < — 1 since some authors argue that K6 is not a distinct spectral type (e.g. 
Kirkpatrick et al.l 119911 ). We follow the convention of omitting K6 as a spectral type and assign 
-1 to K7 and -2 to K5. Thus, we give the spectral types of HIP 29964 and BD -17°6128 A as 
K 5.4 and K5. 7 , resp e ctively, correspond i ng to —1.6 and —1.3 as measured from the relation given 



m 



Reid et al.l (119951 ). IZuckerman et al.l 



2001 



) give the spectral types of these stars as K6/7 and 
K7/M0, respectively, and iTorres et al.l (|2006l ) find K4 and K6. The spectral types o f the stars 
later than K7 agree within the errors with the types given inlZuckerman Song) (120041 ). with the 
exception of HD 155555 C, for which Zuckerman Sz Song ( 20041 ) give M4.5. However, our valu e 
of M3.5 for this star is in agreement with the M3 spectral type found by [Torres et al.l (|2006l ). 
BD — 17°6128 B was not listed independentlv in lZuckerman SongI (j2004l ). but an estimate of its 
spectral type is given in iNeuhauser et al.l (j2002l ) as MO-2, which is close to the spectral type of 
M2.9 that we find here. 



2.2. Equivalent Width of Li I 

Using the splot routine in IRAF, we measured the equivalent widths of all identifiable Li I 
lines at 6707.8A (EW(Li 1) No attempt was made to deblend the Fe 1 hue at 6707.4A from the 
Li I line. This should not affect our later analyses of Li abundances since curves of growth for the 
coolest stars are often derived from measuring the blended lines in model spectra, and thus they 
also include in the quoted EW(Li) values the contribution of the Fe 1 line. The uncertainties given 
reflect the difference between choosing the height of the continuum to be the maxima or the mean 
of nearby features. The Li 1 line is not detected in HD 155555 C, GJ 799 N, or GJ 799 S. Figured] 
shows a comparison of the spectra of the M dwarfs with spectra of main-sequence stars of similar 



■^The spectra of M dwarfs are so rich with hnes that the observed pseudo-continuum is significantly below the true 
continuum level. Because of this, for these stars we actually measure a pseudo-equivalent width rather than a true 
equivalent width. For simplicity, we use the term 'equivalent width' to refer to both equivalent widths (as measured 
for K and earlier spectral types) and pseudo-equivalent widths (as measured for the M dwarfs). 
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Fig. 1. — Spectra of (3 Pictoris Moving Group members in the vicinity of the Li hne (A = 6707.8 
A); the hne's position is marked as a vertical Une. The sohd curve shows the observed spectrum. 
For the M stars, the dotted Une (red in the on-Une edition) is the spectrum of a main-sequence 
template of similar spectral type that has been artificially broadened and shifted to produce the 
best match with the observed spectrum. Filled squares show the difference between the broadened 
template and the observed spectrum, on the same scale as the observations but with an offset of 
0.6 added for clarity. The Li line is clearly detected in the three stars with spectral types earlier 
than M3 (middle row), but it is not present in the later- type stars (bottom row). 
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spectral type that have been rotationally broadened to match the observed projected rotational 
velocity (usini). The difference spectra show that while the line is weak in the early-M stars, it is 
definitely detected, whereas it is clearly not present in the mid-M dwarfs. We measure upper limits 
in the EW(Li) of 30 ml for HD 155555 C, 10 mA for GJ 799 N, and 10 mA for GJ 799 S. 



3. Analysis 



3.1. Luminosity and Rotational Velocity 

We determined the luminosity and projected rotational velocity {v sin i) of each star based 
on data fr om the literature. We c alculated the luminosity of each star using V magnitudes from 
Table 1 in IZuckerman et al.l (j200l|). We assigned a fixed value of cry = ±0.1 as the error on the 
V m agnitudes to account f o r mea surement errors and variability. We used bolometric corrections 
from iKenvon &: HartmannI (jl995l ) appropriate to the spect ral type of the star . Distances to each 
star came from the new reduction of the Hipparcos catalog (|van Leeuwenll2007l ). The V magnitude 
o f BP — 17°6128 is the magnitude of the combined system, so we adopt the luminosities given in 



Neuhauser et al.l (|2002l ) for the A and B components. 



moi\). 



m 



Torres et al 



Torres et al 



Fo r the v sin i o f each star ^ we used the aver age of the values given in IZuckerman et al 
i200d l. and ischolz etaP tooj). For CD 
(|2006l ). vsini = llOkms"^, as 



64° 1 208, w e used the value given 
Zuckerman et al.l (|200lh simply give vsini > 
100 km s-^ BD -17°6128 B does not have a V sin i given in the literature. We measured the 
vsini by artificially broadening a template star spectrum of similar spectral type (Fig. [1]). We find 
that the vsini of BD -17°6128 B is 20 ± 5 km/s. 



3.2. Effective Temperature 



We converted the measured spectral type for each star into an effective temperature Tefj using 
empirical spectral type-Tefj relations to calculate the effective temper ature of each star. For spectral 
types M1-M9, we used the temperature scale given in lLuhmanI (j 19991 ) for "intermediate" stars, since 



our young stars can be exp ected to fall in this class. We combined this with the table given in 
Kenvon &i HartmannI ( 19951 ) for earlier spectral types and interpolated to get effective temperatures 



for our stars. Changes within the spectral type uncertainties lead to typical changes in effective 
temperature of ± ~ lOOK. These data are given in Tabled) 



3.3. Lithium Abundances 



We calculated lithium abundanc es (AjLi) = 12 + I o r( N (L i) /N( }1)) :\ J eSiied 
using the curves of growth given in ISoderblom et al.l (Il993bl ) and IPalla et al. 



2000[) for our stars 
(I2OO7I ) to convert 



Table 1. Late- Type Members of the (3 Pictoris Moving Group 





Spectral 


Distance 


V 


V sin i 


Teff 




EW(Li I) 


A(Li) 




Name 


Other Name Type 


(pc) 


(mag) 


(km s-i) 


(K) 


(Le) 


(mA) 




Notes 



HD 155555 A V824 Ara G5 IV 

HD 155555 B V824 Ara KO IV/V 

HIP 29964 AO Men K5.4 

BD -17°6128 A HD 358623 A K5.7 

CD -64° 1208 •■• K7.2 

HIP 23309 CD -57° 1054 K7.8 

AU Mic CJ 803 M0.7 

GJ 3305 • • • Ml.O 

BD -17°6128 B HD 358623 B M2.9 

HD 155555 C • • • M3.5 

GJ 799 N AT Mic M4.4 

GJ 799 S AT Mic M4.4 



31.4 ±0.5 


7 


21 


34 


57701?% 


31.4 ±0.5 


8 


08 


33 


52501^^^ 


38.6 ±1.3 


9 


77 


15 


4250li« 


45.7 ± 1.6^ 


10 


6"^ 


14 


4130lil° 


28.6 ±0.2 


9 


54 


110 




26.8 ±0.8 


10 


01 


8 


38901J* 


9.9 ±0.1 


8 


81 


9 


37501™ 


29.4 ±0.3 


10 


59 


5 


37101™ 


45.7 ±1.6^ 






20^^ 


3390I70 


31.4 ±0.5° 


12 


71 


6 


33401 


10.7 ±0.4 


11 


02 


10 


321017^! 


10.7 ±0.4 


11 


02 


16 


321017^! 



1 


1 


± 





1 


120 


± 


20 




1,2,4 





54 


± 





06 


250 


± 


60 




1,2,4 





27 


± 





04 


370 


± 


40 


2.5 ±0.3 


1,2,3 





25 


± 





025° 


410 


± 


20 


2.6 ±0.2 


1,2,5 





21 


± 





04 


460 


± 


40 


2.7 ±0.2 


1,2 





15 


± 





03 


370 


± 


50 


2.1 ±0.3 


1,2,3 





07 


± 





01 


120 


± 


60 


-0.1 ±0.7 


1,2,3 





13 


± 





02 


110 


± 


60 


-0.4 ±0.7 


3 





05 


± 





01 


130 


± 


60 


-0.2 ±0.7 


5 





05 


± 





01 




< 


30 


< -1.4 


2 





04 


± 





01 




< 


10 


< -1.4 


2,3 





04 


± 





02 




< 


10 


< -1.4 


2,3 



Not e. — Except where noted, distances are from the updated Hipparcos catalog (jvan Leeuwenll2007^ ■ V magnitudes are from IZuckerman et al 



( 200ll ). and th e other data are o ur ow n measurements. T he given vsini is t he average of the lit erature values cited in the last column. (1) 



given 



Zuckerman et al. ( 2001 ). (2) ^l sin^ given in Torres et al. ( 20061 ) . (3) vsini given in Scholz et al. (2007), (4) Spectral Type from 



Zuckerman et al.l (|2001h . (5) Luminosity from iNeuhauser et al. I (|2002t ) 



''Assuming that BD —17° 6128 A and B have the same parallax as their companion HD 199143. 
''Value for the combined system 

°We give the corrected uncertainty in the luminosity of BD -17°6128 A (|Neuhauser Iboogl ). 

''The vsini of BD — 17°6128 B is not given in the literature. The value given here was estimated from a rotationally broadened stellar template 
^Assuming that HD 155555 C is at the same distance as HD 155555 A/B. 
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the m easured equivalent widths to lithium abundances, following the method used bv lSestito et al 
(|2008l ). We assumed logg = 4.5, which is appropriate for late-type stars descending onto the 
main sequence. We linearly interpolated the curves of growth in A(Li)-EW(Li I) space and used 
a linear extrapolation for abundances < 0.0 dex. For the three stars with effective temperature 
3600K < Tgff < 4000K, we linearly interpolated between the two tables. For the stars for which 
we have only an upper limit to the equivalent width, the calculated abundance is necessarily 
also an upper limit. The quoted uncertainties are the result of individually perturbing both the 
effective temperatures and EW(Li l)s within their uncertainties and adding the effects on the 
derived abundances in quadrature. The results are given in Table [T] 

Because HD 155555 A/B is a spectroscopic binary, we were not able to determine its lithium 
abundance from the EW(Li l)s. The height of the continuum is roughly twice what would be 
measured for each star individually; therefore, the lithium equivalent widths of the individual stars 
would be larger than the equivalent widths g iven in Table [Tj which are measured from the combined 
spectrum. We use the abundances given in lRandich et al.l (|1993) of 3.5 for HP 15555 5 A and 3.7 



for HD 155555 B. Based on the discussion in that paper and in iPasauini et al.l ( 11991 
an uncertainty of ±0.3 dex on those abundances. 



we assume 



3.4. Lithium Depletion Boundary 



We did not detect the LDB of the BPMG. The latest-type stars in our sample, HD 155555 C, 
GJ 799 N, and GJ 799 S, showed no observable lithium lines with very sensitive upper limits (Fig. 
[T]). The LDB of the BPMG has been de tected in two bina ry systems that are members of this group: 
HIP 112312 A/B and V343 Nor A/B dSong et al.ll20o3 . 12003). The 6707.8A Li line was detected 



in the lower-mass stars in these binaries whose spectral types are M4.5 and M5, respectively. The 
stars in our sample with non-detections of the lithium line have slightly earlier spectral types, which 
is consistent with the detected LDB. Comparison of the data with evolutionary models (see below , 
Fig. [3]) implies an LDB age of > 30 Myr. Our data support the conclusion of ISong et al.l (j2002l ) 
that there is a discrepancy betw een the LDB age of th is group and the proposed age from HR 
diagram isochrones of 12~^^ Myr ( Zuckerman et al.ll200ll ). 



3.5. Comparison to Other Work 



To date, two other groups have published lithium abundances and effective temperatures for 
some members of the BPMG (ITorres et al.l l2006l : iMentuch et al.l l2008l ). Given the variation in 
effective temperatures and lithium abundances measured by each group, our data are consistent 
with previous measurements. Overall, we see the same trend as other groups with lithium decreasing 
as a function of decreasing spectral type and only upper limits in lithium abundance measured for 
GJ 799 N and S. 
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3.6. Comments on Individual Objects 



3.6.1. CD -64° 1208 



Examination of the 6707.8A Li I line and the Ca I hne at 6717.7A suggests that CD -64° 1208 is 
a double-hned spectroscopic binary (Fig. [2|) . Cross-correlation of the spectrum with radial velocity 
standards of known velocity and similar spectral type gives velocities of — 28 ± 8 km/s for the 
secondary and 43 it 5 km/s for the primary. Based on the height of the cross-corre lation peaks, 
the se condary has roughly 70-80% the luminosity of the primary. We note that [Torres et al 



(1200a) also flag this system as "SB?", i.e. a po ssible spectroscopic binary. Our derived spectral 
type of K7 agrees with that of iRiaz et al.l (|200g) and is similar to the K5 spectral type found by 
Torres et al.l (|2006l ). We note that the vsini of 110 km/s measured bv I Tor res et al.l (|2006l ) is likely 
an overestimate. The rotation is rapid, but some of the line broadening is due to the velocity 
separation of the two stars. If the luminosities and effective temperatures of the two components 
are similar, as seems likely from the spectrum, the equivalent widths of the individual lithium line 
and the lithium abundances of the individual components should be roughly the same as those 
given in Tabled) 



3.6.2. AU Mic 

If AU Mic is the same age as the other members of the BPMG, it appears to be under-luminous 
in the HR diagram (Fig. [3]). We observe log(L/L0) = —1.15, but from the figure we might expect 
a value closer to log(-L/L0) = —0.95 for it to lie on the same HR diagr am isochrone as th e other 



BPMG members. AU Mic is known to have an edge-on debris disk (cf. iKalas et al.ll2004l ). which 
could lead to some extinction. Given the difference between our observed and expected values for 
the luminosity, if AU Mic is coeval on the HR diagram, then we estimate an extinction due to the 
debris disk of Ay ~ 0.5 mag. Since this extinction is not well characterized, and it is dependent on 
an assumption of coevality, we do not take it into account in our analysis. This proposed extinction 
should not affect our measurement of the effective temperature since that is determined from the 
spectral type which is measured over a small range in wavelength. However, it does affect any age 
determined from the HR diagram position. If the luminosity we observe is less than the intrinsic 
luminosity of AU Mic, we would overestimate the age of the star (see Fig. [3]). 



4. Comparison to Evolutionary Models 

In Figu r e [3] we plot our data agains t three sets of evo lutiona ry models: D'Antona &: Mazzitelh 



(119971 . Il998|), iBaraffe et all (|l998l ). and ISiess et all (|2000l ). The ISiess et al.l tOOd ) models include 
the effects of convective overshooting. Known multiple systems are plotted in color. We also plot 
HIP 112312 A/B and V343 Nor A/B, two systems that cross the LDB. We calculated luminosities 
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Fig. 2.— Spectrum of CD -64° 1208 showing the Li 1 (6707.8 A) and Ca I (6717.7 A) hnes. Note 
that each hne appears to have two components, shown by tick marks at measured RVs of —28 and 
43 km/s. 
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Fig. 3. — Comparison of the data and theore t ical models. The d ata ar e overlaid on thre e sets 
of models: IP'Antona fc Mazzitellil (|l997l . Il998l ). Isaraffe et alJ (|l998l ). and Isiess et aP hood ). The 
results of this work are shown as diamonds. Single stars are shown in black and systems with 
multiple stars are shown in color (HD 155555 A/B/C (green); BD -17°6128 A/B (orange); and AU 
Mic / GJ 799 N/ S (red)). Upper l imits are shown as open diamonds with arrows. The square s show 
HD 112312 A/B dSong et al.ll2002h and the triangles show V343 Nor A/B dTorres et ahlbood ^. Note 
that the luminosities and effective temperatures for GJ 799 N/S are similar, so the two stars appear 
to be a single data point. Likewise, the lithium abundances and temperatures for HD 112312 B 
and V343 Nor B are the same. Notice that the ages predicted by the location on the HR diagram 
and by the location on the corresponding lithium abundance plot differ for any individual star. 
Taking the stars as a group, it is apparent that the models do not correctly predict the observed 
patterns of lithium depletion. This is also true of the stars in individual multiple systems whose 
components must be coeval. 
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and effective temperatures for th ese stars in the same rnanner as for our sample, using pub lished 



spectral types and V magnitudes (IZuckerman et al.ll200ll : ISong et al.ll2003l: [Torres et al 



used the lithium abundances from Torres et al 



((20061). Neither ISong et alJ ((2002) nor lTorres et al 



2006). We 



(j2006l ) quote a lithium abundance for HIP 112312 A. S ince t he spectral type and the upper limit 
in the equivalent width of this star given in lSong et al.l (|2002l ) are similar to those of GJ 799 N, we 
assign HIP 1123112 A an upper limit of A(Li) = —1.4 dex. 

To explore the differences between the models and our data, we independently determine ages 
for each star in two ways: from its position on the HR diagram and from its position relative to the 
lithium depletion isochrones. We estimate the uncertainty of the HR diagram age by perturbing the 
effective temperatures and luminosities within their uncertainties and adding the resulting age errors 
in quadrature. Because the measured lithium abundance is correlated with effective temperature, 
we estimate the uncertainties in the lithium depletion age by simultaneously varying the effective 
temperatures and lithium abundances within their uncertainties and taking the extreme values for 
the age as a measure of its uncertainty. These data are given in Tableland shown in Figure HI This 
figure shows that the ages calculated from the HR diagram are roughly consistent with a single age, 
but the ages measured from the lithiu m depletion have systernatic trends. The trend of inc r easing 
age with spectral type found using the iD'Antona &: Mazzitellil (jl997l . Il998l ) and ISiess et al.l (j2000l ) 
lithium depletion isochrones simply reflects the fact that t he data cross t he iso chrones instead 
of lying parallel to them, as seen in Figure [3l Although the iBaraffe et al.l (|l998l ) models are the 
closest to predicting a coeval group using lithium depletion, they still fail to accurately represent the 
latest-type stars in th e sample. Even thou gh the predicted ages are similar, the pattern of lithium 
depletion seen in the iBarafFe et al.l (|l998l ) models is the opposite of what is seen in our data. In 
particular, these lithium depletion isochrones indicate that the lithium depletion boundary occurs 
between 3400K and 3200K for the 30 Myr isochrone, which approximates our data. Thus, we 
expect that stars at the cooler end of the temperature range 3200-3400K should have a higher 
lithium abundance than stars at the warmer end. In contrast, we observe that HD 155555 C, GJ 
799 N, and GJ 799 S with rcff~3300K, 3200K, and 3200K, respectively, have only upper hmits in 
lithium abundance, whereas the warmer star BD — 17°6128 B (Tcff~ 3400K) has a larger lithium 
abundance. 



From these two methods for measuring ages, we observe two systematic trends in the data. 
First, the ages of the mid-M dwarfs calculated from the HR diagram isochr ones are slightly youn ger 
than the ages of the warmer stars. This effect has also been noted by IStassun et al.l (|2004l ) in 
the ages of binary com ponents, and has been r ecognized more generally as a problem for PMS 
evolutionary models by [Hillenbrand et al.l (|2008l ). The second trend is that the ages of the mid-M 
stars calculated from their lithium depletion are older than the ages of the warmer stars. 

While the HR diagram isochrones match our data well and the stars are roughly coeval, none 
of the evolutionary models have lithium depletion isochrones consistent with a coeval group. If 
we consider each multiple star system independently of the group as a whole, we find similar 
inconsistencies in the predicted ages for the hotter and cooler stars in each binary. Since the two 
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stars in each binary are almost certain to be coeval within < 1 this indicates that the 

mismatch between our data and the isochrones is not the result of contamination by unrelated 
stars or an extended period of star formation. In general, the models seem to have difficulty 
describing the lithium depletion of the lowest-mass stars. In particular, the modeled rise in lithium 
abundance at lower temperatures (the LDB) occurs at temperatures too high for our stars in all 
of the evolutionary models. Given the nominal HR diagram age, this means that an estimate of 
the age of the group from the LDB will be systematically older than an age from fitting the HR 
diagram regardless of the model used. 



5. Larger M-dwarf radii to the rescue? 

We have shown that the pattern of Li depletion seen in the BPMG stars is at odds with the 
predictions of PMS Li depletion models, in the sense that the models under-predict the amount 
of Li depletion in the M stars in our sample. The evolutionary models have a certain amount 
of intrinsic uncertainty due to assumptions about mixing length, equation of state, opacity, and 
atmosph ere and boundary conditions, whose effects on lithium depletion are discussed in some 



detail in iBurke et al.l (120041 ). Without a specific physical reason to change the assumptions, it is 
difficult to decide which parameter or combination of parameters should be changed in order to 
fit the lithium data. However, there is an increasingly large body of observations (see below) that 
evolutionary models tend to under-predict the radii of low-mass stars, particularly M stars, which 
may point the way toward the types of modifications to the models necessary to bring the HR 
diagram ages and Li depletion ages into agreement. 

Changes in stellar radii could affect the observed pattern of Li depletion in two ways: by 
actually changing Li depletion as a function of mass, and/or by altering the effective temperatures 
of stars, thus shifting them relat i ve to the model isochrones. Regarding actual changes in Li 
depletion, the work of King et al. (j2009l ) suggests a connection between Li depletion and stellar 



radii. They propose that a range in radii at a given stellar mass (perhaps due to a range in rotation 
rates and/or chromospheric activity) would lead to a range of interior temperature profiles and 
could explain the dispersion in lithium depletion in the Pleiades. 

However, even if the models are essentially correct in terms of the interior conditions for a given 
mass, which largely govern the rate of Li depletion, there could still be a disagreement between the 
model predictions and the observations if the models do not correctly predict the surface conditions, 
e.g. the radius or effective temperature. We suggest here that current data show evidence of exactly 
this sort of disagreement. 



^The multiplicity frequen cy of stars at the you ngest observable pre-main-sequence ages is comparable to or greater 
than that among field stars (jPuchene et al.ll2007l . and references therein), indicating that binary formation is part of 
the star formation process. Thus, we can consider binary systems to b e N = 2, coeval clusters to within < 1 Myr, 



although there are hints of smaller age differences ( Stassun et al. 20081 ) 
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Fig. 4. — Age as a function of spectral type, calculated for individual stars from different models. 
The left column shows ages calculated from the HR diagram and the right column shows ages 
calculated from the lithium depletion isochrones. Since GJ 799 N and S have the same spectral 
type, the points for GJ 799 S have been displaced by +0.1 in spectral type for clarity. 
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Fig. 5. — iBaraffe et alJ (|l998l ) evolutionary models with artificially inflated radii. The solid lines 



show the isochrones with an ad hoc correction assuming that for Tcff<4200K the radii should be 
increased by 10%. The dotted lines show the original model isochrones from Fig. [3j The other 
symbols are the same as those shown in Fig. [3l 
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There is mounting observational evidence that main-sequence stellar evolutionary models pre- 
dict radii for M dwarfs that are 10-20% smaller than the observed radii. This effect is seen 



in stellar angular diameters measured interferometrically (I Lane et al 



2008 



Berger et alj 120061) . in radii determined from eclipsing binaries (jLopez-Morales ,2007i : iRibas et al 



2001 



Segransan et al. 



2003 



Fernandez et al.ll2009l ). and i n radii determined from detailed modeling of multi- band pho 



tometry to measure LboI and Teg (jMullan h MacDonald 



2001 



Casagrande et alJ 120081 ) . There is 



still some qu estion as to whether this discrepancy exists only for s tars that are the most m agneti- 
cally active ([Pemorv et al.ll2009l ). or if it applies to all M dwarfs (ICasagrande et al.ll2008l). There 



does seem to be at least some differential effect from stellar activity (jMullan &: MacDonaldl 12001 



Morales et al.ll2008l ). and m odels incorporating t he effects of stellar activity are better able to re- 
produce the observed radii (jChabrier et alJ 120071 ). For the present work, the question is moot as 
the stars in our sample (and indeed all late-type pre-main-sequence stars) are very active, with 
large X-ray luminosities and active chromospheres. 

The discrepancy between the true radii and the radii predicted by models very likely exists 
in the pre-main-sequence phase as well, if it is due to effects of magnetic activity and/or missing 
opacity sources in the models. If this discrepancy is present, how would it manifest itself? Here 
we argue that the effect would be similar to what we observe, and correcting for this would bring 
HR diagram ages and Li depletion ages closer together. Corrected model isochrones for a given age 
would need to move to lower Tgff , increasing the age inferred from the HR diagram and decreasing 
the age inferred from Li depletion for a given star. 

First consider the HR diagram. Pre-main-sequence stars change position in the HR diagram 
as they age precisely because their radii are changing, as the stars descend (contract) toward the 
main sequence. If, at any given age, the true radii of stars are larger, then the model isochrones 
for that age need to move up and to the right. Thus, the HR diagram age of a given star (of 
fixed, observed Tgfj and -Lroi) inferred from these corrected models would be older than before. 
Put another way, if radii are inflated by effects other than youth, it will take longer for a star to 
contract to a given radius, and thus to reach a given position on the HR diagram. 

Now, consider the A{L\)-T(,s plane. The amount of Li depletion at a given age is primarily 
driven by a star's mass, which sets its interior conditions. Inflating the radius of a star will cause 
that star to have a cooler surface, i.e. a lower T^s- Thus, corrected model isochrones would need 
to shift to the right for a given j4(Li) and age. Put another way, a star with an inflated radius will 
have a hotter interior than we might otherwise expect given its surface appearance, and thus will 
exhibit more Li depletion at a given age and T^^ than would be predicted if the radius were not 
modified. This matches well with the enhanced Li depletion we observe for the cooler stars. 

Calculating new pre-main-sequence evolutionary models incorporating these effects is beyond 
the scope of this work, but we can explore these effects to a limited degree by introducing an ad 
ho c shift to the e xistin g model isochrones. As a test of this idea, we shifted the model isochrones 
of iBaraffe et al.l (| 19981 ) and compared them to our observations of A(Li) and Tgfj for the BPMG 
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stars. Based on the observations of M stars di scussed above, we assum ed that the larger radii are 
present only for stars cooler than ~ 4200 K (jCasagrande et al.l 120081 ) . and that the model radii 
should be inflated by 10%. As shown in Figure O these corrections to the models do indeed bring 
the inferred Li depletion and HR diagram ages of the cooler stars into better agreement, though 
some discrepancies remain. Clearly, a more systematic treatment of this effect is needed, perhaps 
taking into account effects of different activity levels as well as different masses. 

To close, we note that one implication of this result is that neither the LDB ages nor the HR 
diagram ages of M stars are correct, but that the answer lies somewhere between the two. 



6. Conclusions 



We have measured lithium abundances and effective temperatures for ten members of the 
(3 Pictoris Moving Group. We compare these abundances to the p redicti ons of the e voluti onary 
models of IP'Antona fc Mazzitellil (Il997l . \l9m), Isaraffe et aP dlQQsl ^. and Isiess etliD (|20od ). We 
find that while the HR diagram isochrones of these models reproduce our data fairly well, the 
lithium depletion isochrones do not. In particular, the models predict less lithium depletion for the 
latest-type stars than we observe. Thus, the ages determined based on the lithium abundances are 
older than the ages determined from the HR diagram. 

The most striking result is the non-detection of lithium in the latest-type members (M3.5-M5) 
of this group. The lack of an observable lithium line in HD 155555 C, GJ 799 N, and GJ 799 S 
indicates that lithium depletion proceeds at a much faster rate than is predicted in the theoretical 
models given their spectral types. One potential clue lies in the discrepancy between the observed 
radii of M dwarfs and the radii predicted by models. 

Regardless of the physical mechanism, it appears that mid-M stars deplete their lithium rapidly. 
The presence of the lithium line at 6707. SA in a star of this type remains a clear indicator of youth 
and suggests an age younger than that of the /3 Pictoris Moving Group, while the lack of a detectable 
line is not a strong indicator that a star is old. 
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Table 2. Ages in Myr Calculated from Different Models and Isochrones 
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